clear
clc
N1 = 600
N2 = 75
L = 0.01.*47.123
M = 0.01
S = 0.0001.*1.3273
u_r = 4.*pi.*power(10,-7)
load data.mat

table2 = zeros(3,6)
table2 = M2
H_0 = (2.*sqrt(2).*N1.*M2(1,:))./(L.*M2(2,:))
table2(3,:) = H_0
table3 = zeros(3,6)
table3= M3
B_0 = (2.*sqrt(2).*M.*M3(1,:))./(N2.*S.*M3(2,:))
table3(3,:) = B_0
AverageH_0 = mean(H_0)
AverageB_0 = mean(B_0)

table1 = zeros(6,10)
table1 = M1
table1(4,:) = AverageH_0 .* M1(2,:)
table1(5,:) = AverageB_0 .* M1(3,:)
table1(6,:) = table1(5,:) ./ (table1(4,:).*u_r)

T(1,:) = table1(1,:)
T(2,:) = table1(6,:)
[n, m] = size(T)
x = T(1,:);
y = T(2,:);
plot(x,y,'o');
figure(1)
U = 0
u_r0 = -0.5557.*U.^2 + 61.72.*U + 509.3